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FOREWORD 


This  report  was  prepared  for  Applied  Mathematics  Research  Laboratory, 
Aerospace  Research  Laboratories  by  D.  S.  Clemm,  P.  R.  Krishnaiah  and  V.  B. 
Waikar  under  Project  7071,  "Research  in  Applied  Mathematics".  The  work  of 
V.  B.  Waikar  was  performed  at  the  Aerospace  Research  Laboratories  while  in 
the  capacity  of  an  Ohio  State  University  Research  Foundation  Visiting 
Research  Associate  under  Contract  F  33615-67-C-1758.  The  present  affiliation 
of  V.  B.  Waikar  is  Miami  University,  Oxford,  Ohio. 

In  this  report,  the  authors  gave  tables  for  the  exact  percentage  points 
of  the  joint  distribution  of  the  extreme  roots  as  well  as  the  tables  for  the 
exact  lower  percentage  points  of  the  largest  root  of  the  Wishart  matrix. 

The  authors  wish  to  thank  Miss  Eva  Brandenburg  for  typing  the  manuscript. 


ABSTRACT 


Let  and  2  be  respectively  the  smallest  and  largest  roots  of  the 
central  pxp  Wishart  matrix  with  n  degrees  of  freedom.  In  this  report,  the 
authors  gave  tables  for  the  exact  values  of  U  for  p  *  2(1)10(2)20,  a  =  0.05, 
0.025,  0.01,  0.005  and  different  values  of  n  where 

P[U_1  $  U]  =  (l-2a). 

Also,  exact  lower  II,  2.51,  SI,  101  points  of  the  distribution  of  2p  are 
given  for  p  =  2(1)10(2)20  and  for  different  values  of  n. 
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I.  INTRODUCTION 


Krishnaiah  and  Chang  [3]  gave  an  expression  for  the  probability  inte¬ 
gral  associated  with  the  joint  density  of  the  extreme  roots  of  the  Wishart 
matrix.  For  computational  purposes,  this  expression  is  significantly 
better  than  other  expressions  given  in  the  literature  Ql],  16]).  Using 
this  expression,  the  authors  constructed  tables  for  the  percentage  points 
associated  with  the  joint  density  of  the  extreme  roots  of  the  Wishart 
matrix  by  imposing  the  condition  that  the  lower  critical  value  is  equal  to 
the  reciprocal  of  the  upper  critical  value.  These  tables  are  useful  in 
the  application  of  Roy's  test  for  testing  the  hypothesis  that  the  covariance 
matrix  of  a  multivariate  normal  population  is  equal  to  a  specified  value. 
Some  computations  are  also  made  to  check  for  the  accuracy  of  the  approximate 
tables  constructed  by  Hanumara  and  Thompson  [2]  by  putting  a  different 
restriction  on  the  choice  of  the  critical  values.  Also  tables  of  the  lower 
percentage  points  of  the  largest  root  of  the  Wishart  matrix  are  constructed 
in  this  paper. 

2.  PROBABILITY  INTEGRAL  OF  THE  JOINT  DENSITY  OF  TIE  EXTREME  ROOTS 

Let  X  be  a  pxn  matrix  whose  columns  are  independently  distributed  as 

multivariate  normal  with  zero  mean  vector  and  covariance  matrix  I  ,the 

P 

identity  matrix.  Also  let  S  *  XX'  and  let  *■•••*■  £p  be  the  latent 

roots  of  S.  Then  it  is  well  known  that  S  has  Wishart  distribution  and  the 
joint  density  of  £j,  ...,  £p  is  given  by 


(2.1) 


r 


P  r  p 

fiUi.-.-.O  =  k(p,n)  n  (i,  exp(-*./2)  n  (*■-  i.) 

P  i»l  1  x  i>j  1  J 


where  r  =  (n-p-l)/2  and 


k(p,n)  =  »p/2(l/2)nP/2  n  [r(n+l-i)/2) r((p+l-i)/2] . 

i=l 

Krishnaiah  and  Chang  [3]  gave  an  exact  expression  for  the  probability  integral 
of  the  joint  density  of  and  4.^  which  is  given  by 

P[L  .<  <  U]  =  k(p,n)o(f;  p,r ,L,U)  (2.2) 


where 


and 


r 

aCy;  2m,r,L,U)  when  p  =  2m 
oU;  p.r.L.U)  2m 

■  t  (‘ljl  Fr*iLW  2m+l,r ,L,U) 


(2.3) 


i=0 


i+] 


when  p  «  2m+l 


1/2 


(2.4) 


A(y;  2m,r,L,U)  *  ((a^)^ 

Gt(T;  2m+  1,r,L,U)  =  1  (aij)i, j*l . t-1,  t+1 . 2*1 

aij  -  fs  -  Fs  •  1  . 


1/2 

i  for  p>  1 


Fs  =  5  Fs(e)  0t  'i'{0)de  *  Fs(0)  *  J6  xS  to  » 

’i'(x)  =  exp(-x/2).  Also,  let  Gjfo  l,r,L,U)=l. 

In  this  paper,  the  authors  give  tables  for  the  exact  values  of  U  for 
«■  0*05,  0.025,0.01,  0.005,  p  =  2(1)10(2)20,  n  *  (p*l) (1)20(2) 30(5) 50  where 

P[U_1  £lx  i  ip  i  U)  =  (l-2a).  (2.5) 


2 


The  authors  also  gave  lorer  1%,  2.5$,  5$  and  10%  values  of  the  distri¬ 


bution  of  i  for  p  -  2(1)10(2)20,  n  =  (p+1)  (1)20(2)30(5)50  by  using  the 
following  known  formula  [3]: 

P[£p  >  U]  «  l-k(p,n)P(u<;  p,r,0,U) 

llanunara  and  Thompson  [2]  constructed  tables  for  L  and  U  for  p  =  2(1) 

10  and  various  values  of  r  and  a  satisfying 

P[  L  i  <  tp  <  U]  =  l-2a  (2.6) 

and 

P[  tx  >.  L]  =  1  -  a.  (2.7) 

For  p  =  2,  Thompson  [7]  computed  exact  values  of  L  and  U.  For  p>3, 

Hanunara  and  Thompson  [2]  state  that  exact  values  are  very  difficult  to 
compute  and  so  they  computed  approximate  values  and  Uj  of  L  and  U 
respectively  where  Lj  and  satisfy 
V[  l.  »  LJ  *  ri  -  0) 

1  1  (2.8) 

P(  <  Uj]  =  (1  -  a). 

For  p  -  2  and  3,  they  found  empirically  that  and  Uj  are  close  to  L  and  U 
respectively.  For  higher  values  of  pfg  10)  we  computed  L  and  U  exactly  for 
some  typical  values  of  n  and  a  and  found  them  to  be  close  to  and  given 
in  [2] . 


3.  APFLICATKNS 

Let  IL;  I  =  E„,  H  :  E  #  E_  where  E  is  the  covariance  matrix  of  a 
o  0’  a  u 

p-variate  normal  distribution  and  E^  is  kriown.  Further  let  S/n  be  the  sample 
covariance  matrix  based  on  a  sample  of  size  n+1.  According  to  the  procedure  of 


Roy  [5,  p.  30]  an  acceptance  region  for  testing  MQ  against  II  is  given  by 


L  *  *  U 


(3.1) 


-1 


where  i .  <  *,<  ...  <  t  are  the  characteristic  roots  of  S£„  and  L<U  are 

1  Z  p  vJ 

constants  such  that 


P(  L  {  <  I  i  Uj  lip]  =  1  -  2a. 


(3.2) 


Note  that  under  Hq,  SEq1  has  a  Wishart  distribution  W(I,n).  Now,  the  optimum 
choice  of  L  and  U  (in  the  sense  of  maximizing  the  power)  is  not  known.  Thus 
one  simple  way  of  choosing  L  and  U  is  to  put  L  =  1/U  and  then  use  the  tables 
given  in  this  paper. 

Che  can  also  construct  simultaneous  confidence  intervals  for  the  elements 
CTj  of  the  covariance  matrix  r  using  the  tables  given  here.  Roy  [5,  p.  106] 
showed  that 


if1  a,'  S  a,  *  a,'  IaN<  L'1  a'  Sa 


(3.3) 


is  a  set  of  simultaneous  confidence  bounds  on  a'  Ea,  for  all  arbitrary 
nonnull  vectors  a  with  confidence  coefficients  1  -  2a  provided  L  and  U  satisfy 
equation  (2.6).  By  choosing  a^  appropriately  in  (3.3)  one  can  obtain  simul¬ 
taneous  confidence  intervals  for  Thompscn  [7]  derived  the  following 

simultaneous  confidence  intervals  for  cm  .  using  Roy's  results: 


If1  s..*  oa.<  L*1  sa  i  =  1.....P 


\°n  -  r  (U'1-  L’^s.J  i  \  (L'1-  if1)  (s..  s.  •) 


1/2 


(3.4) 


13 


i)1 


li 


i*j. 


Thus  (3.4)  represents  simultaneous  confidence  intervals  for  cm..,  i,j*l,...,p 


4 


with  confidence  coefficient  1  -  2a  where  L  and  U  satisfy  (2.6)  and  S 


(s^).  Hence  putting  L  =  1/U  and  using  Table  I,  one  can  construct  exact 
simultaneous  confidence  intervals  for  given  by  (3.4). 

As  a  further  application,  we  mention  that  these  tables  are  useful  in 
obtaining  the  simultaneous  confidence  intervals  for  the  variance  components 
of  the  two-way  layout  with  unequal  variances.  Since  this  application  is 
discussed  in  detail  by  Thompson  [7,8]  and  Ilanunara  and  Thompson  [2],  we 
omit  its  discussion  here. 


I 


1 
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